Role of defects in the electronic properties of amorphous/crystalline Si interface 
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The mechanism determining the band alignment of the amorphous/crystalline Si heterostructures is 
addressed with direct atomistic simulations of the interface performed using a hierarchical combina- 
tion of various computational schemes ranging from classical model-potential molecular dynamics to 
ab-initio methods. We found that in coordination defect-free samples the band alignment is almost 
vanishing and independent on interface details. In defect-rich samples, instead, the band alignment 
is sizeably different with respect to the defect-free case, but, remarkably, almost independent on the 
concentration of defects. We rationalize these findings within the theory of semiconductor interfaces. 
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Amorphous-crystalline silicon interfaces (a-Si/c-Si and 
a-Si:H/c-Si) are of unquestionable technological impor- 
tance for their use in solar cells and in other optoelec- 
tronic devices. Their electronic and optical properties 
are governed by the band discontinuities whose 
value and origin however have not been clearly estab- 
lished so far. To our knowledge, few and conflicting ex- 
perimental results exist for a-Si:H/c-Si, ranging from 
to 0.7 eV for the valence band offset (VBO) (see Refs. 

and references therein), and even less informations 
are available for a-Si/c-Si 0. On the theoretical side, 
the problem of band alignment is addressed in two pa- 
pers ]5],[7j , but only from a comparison of the two phases 
separately and not from a direct simulation of the in- 
terface. The atomic-scale modelization of the interface 
and its structural properties are the subject of few other 
works |§|-|rT| which, however, do not examine the elec- 
tronic structure. The combined study of the structural 
and electronic properties is necessary to understand the 
mechanism determining the band alignment and in par- 
ticular the role of bulk-related and interface-specific de- 
tails. 

With this aim we address here our attention to the 
a-Si/c-Si junction, described using periodically repeated 
supercells. To generate the structures we adopted 
a hierarchical combination of models and simulation 
schemes, including the well known Wooten- Winer- Weaire 
(WWW) model Q, model-potential Jl3| and semi- 
empirical tight-binding (TB) molecular dynamics (MD) 
pq |. We used then ab initio state-of-the-art pseudopo- 
tcntial calculations Jl6|,[l7j to refine the equilibrium ge- 
ometries and to compute the electronic structure. 

We will consider a-Si/c-Si structures with different a- 
Si atomic-scale morphologies. In order to investigate 
the bulk-related effects, we will also study for compar- 



* E-mail: peressi@ts.infn.it 



ison the band lineup at crystalline strained/unstrained 
Si homojunctions and at heterocrystalline cubic dia- 
mond/hexagonal diamond interfaces. 

Atomic-scale model for a-Si/c-Si interface. - For the 
description of the (OOl)-oriented a-Si/c-Si interface we 
consider periodically repeated tetragonal (001) supercells 
containing 320 Si atoms, with cell dimensions a = 2oq 
and c ~ 10ao, where ao is the lattice parameter of the 
crystalline phase. The analysis of the electronic struc- 
ture, in particular of the spatial decay in the crystalline 
region of the gap states induced by coordination defects 
of the a-Si region, shows that this cell size is large enough 
to recover crystalline bulk features in the middle of the 
c-Si slab. 

The atomistic models of the interface structure are ob- 
tained according to a procedure that, basically, consists 
in the following four steps: (i) a tetragonal cell of a- 
Si with dimensions a = 2oq, c — 5ao and 160-atoms is 
created by quenching from the melt during a constant- 
pressure MD simulation based on the "environment de- 
pendent interatomic potential" (EDIP) jl3); the quench- 
ing procedure was operated by velocity rescaling at a rate 
as slow as 0.1 K/ps (i.e. about 100 times slower than typi- 
cal TB-MD simulations Jl4| ) , so that density fluctuations 
during the liquid-to-amorphous transition were properly 
damped and a fully stress-free a-Si sample was indeed 
obtained; a sample interface is then created by joining 
through the smallest face the above a-Si sample to a c-Si 
one with the same cell; (ii) the structure is further an- 
nealed up to 800 K (heating rate ~ 1 K/ps), aged at that 
temperature for about 100 ps, and finally relaxed down 
to K (cooling rate ~ 1 K/ps) by means of constant- 
pressure, constant-temperature EDIP-MD simulations; 
(in) at this point we switch from EDIP-MD to TB 
MD, adopting the parametrization described in Ref. Jl8| ; 
since the two models give different equilibrium lattice 
parameter for c-Si (o,q DIP =5.43 A to be compared with 
a TB =5A5 A ), we first perform an homogeneous rescal- 
ing of the supercell dimensions and atomic positions, thus 
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keeping the optimal c/a ratio provided by the previous 
two steps; (iv) the internal degrees of freedom are then 
further relaxed by means of a new constant-volume TB- 
MD annealing procedure (annealing temperature is 1000 
K, the heating/cooling rates are, respectively, 500 K/ps 
and 160 K/ps). Overall the computational workload is 
heavy - especially for the TB-MD steps - and has to 
be mastered by means of O(N) formulation of TB-MD 
and/or extensive use of distributed computing 

The resulting structure is then used as input of ab 
initio calculations, performed within the framework of 
density-functional theory in the local density approxi- 
mation, using norm-conserving pseudopotential |l9f| , a 
plane-wave cutoff of 16 Ry, and 4 k points generated by 
a symmetrized Monkhorst-Pack grid |po| . We first per- 
form a further homogeneous rescaling of the supercell di- 
mensions and atomic positions according to the new the- 
oretical equilibrium lattice parameter ao CF =5.47 A of 
c-Si. We also perform an additional optimization of the 
internal degrees of freedom through minimization of the 
total-energy and atomic forces, and finally we use the re- 
sulting configuration for the calculation of the electronic- 
structure. Because of the non-periodicity in the amor- 
phous slab (see Fig. Q), it is not possible to define an 
average electrostatic potential in that region and there- 
fore to extract the band alignment according to the usual 
procedure explained in Ref. |21| for the crystalline junc- 
tions. Instead, from the supercell calculations we can ex- 
tract the Fermi energy Ep and the average electrostatic 
potential in the crystalline region, (V) c , both measured 
with respect to the average electrostatic potential of the 
entire supercell which is set to zero (see Fig. [2]). From 
the electronic structure calculations of bulk c-Si we ob- 
tain the position of its valence band top edge E VjC with 
respect to its own average electrostatic potential, and 
therefore we can calculate the band alignment in terms 
of the Schottky barrier: <fi p = Ep — E v , c — {V) c . The nu- 
merical uncertainty associated to (f> p , coming mainly from 
the determination of Ep, is of the order of 0.1 eV. Cal- 
culating the band alignment in terms of (f> p is consistent 
also with the fact that we treat the a-Si/c-Si interface as 
a metallic system, due to the presence — in general — of 
defect-induced gap states. 

We discuss the results for 6 different a-Si/c-Si config- 
urations, obtained by means of different details in the 
above preparation procedure. Variations on such a pro- 
tocol were aimed at producing a-Si samples with varying 
concentration of coordination defects. In Tab. I we re- 
port the relative occurrence in the amorphous part of 
under-coordinated and over-coordinated Si atoms, calcu- 
lated simply using a bond-cutoff distance 15% larger than 
the crystalline bond length. This choice is consistent with 
the actual position of the first minimum of the pair cor- 
relation function computed in the amorphous slab. The 
six samples are representative of at least three quite dif- 
ferent atomic-scale arrangements, with very low (sample 



A), intermediate (B, C) and high (D, E, F) concentra- 
tion of defects. Such a diversity of relative occurrence of 
coordination defects also reflects in the absolute average 
atomic coordination that varies from 3.99 for sample A 
(small under-coordination) to 4.3 for sample F (sizeable 
over-coordination). These values are consistent with ex- 
perimental findings on real a-Si samples |Q. Within the 
last two categories (samples B, C and D, E, F, respec- 
tively), we present the results for two or three samples 
in order to examine the possible influence of fine struc- 
tural details. The cell dimension c/a is slightly varying 
in the different samples, with positive and negative small 
deviations from the ideal value c/a=5. 

Results for deject-free a-Si/c-Si samples. - We first 
focus on sample A, which, at variance with our typical 
procedure, has been obtained starting directly from a 
320-atom supercell and not joining two separate a-Si and 
c-Si cells as described in step (z). The amorphization of 
the a-Si slab has been achieved by the WWW model ; 
then the sample has been treated with EDIP-MD for the 
volume optimization, and with a TB athermic annealing 
for a final optimization of the atomic positions. Sam- 
ple A has only two three-fold coordinated defects, one 
at each interface. It has a clear semiconducting charac- 
ter, with a well defined energy gap, and a band align- 
ment 4> p =0.06 eV. In this case it is possible to calculate 
also the VBO, since the position of the topmost valence 
band edge in the amorphous slab with respect to Ep is 
clearly obtained from supercell calculations. We found 
VBO=0.04 eV, with c-Si higher, sizeably smaller than 
the predictions of Refs. §,§ (0.25 eV and 0.19 eV re- 
spectively). The discrepancy may be ascribed to the fact 
that those calculations do not consider directly the inter- 
face; in Ref. |?J, for instance, the electrostatic potential 
lineup is assumed to be zero. 

Because of its constituents, sample A can be re- 
garded as an isovalent lattice-matched semiconduc- 
tor/semiconductor crystalline heterojunction. Indeed it 
obeys the general trend of this class of junctions, where 
the VBO is a bulk- related effect, independent on the in- 
terface details j2^] : remarkably, the two interfaces in the 
supercell, although slightly morphologically different, are 
really equivalent in terms of electronic structure, since no 
electric fields are present in the bulk slabs. 

This behavior is consistent with experimental findings 
for the analogous system a-Si:H/c-Si m, where a constant 
VBO=0.44 eV has been found for different samples char- 
acterized by the same amount of H but different thickness 
of the interface region. 1 22 

Band alignment at Si/Si crystalline junctions. - In 
order to examine the role played by the density and the 
crystalline symmetry on the band alignment, wc consider 
two particular crystalline Si junctions: the unstrained di- 
amond/strained diamond and the hcterocrystalline cubic 
diamond/hexagonal diamond interfaces. For the former 
case we consider specifically a positive strain of 1% along 
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the (111) direction, for a direct comparison with the 
case of hexagonal diamond which naturally exhibits such 
strain: we found VBO=0.05±0.01 eV (with the strained 
slab higher) after optimization of the internal degrees of 
freedom in the strained Si slab. The contribution to the 
VBO due to the change in density is 0.02 eV, related to 
the hydrostatic deformation potential p3[ . The remain- 
ing contribution comes from the strain-induced splitting 
of the valence band top edge (0.10 eV). 

We simulate the 

heterocrystalline cubic diamond/hexagonal diamond Si 
junction with a supercell with 30 atoms, namely with 
3 unit slabs for both phases along the (111) direction 
|pif . After a careful optimization of the macroscopic and 
microscopic degrees of freedom, which gives for the struc- 
ture of the hexagonal diamond results close to those of 
Rcf. |H, we found VBO=0.23±0.01 eV, with hexagonal 
diamond higher than cubic diamond, as in Ref. j26|. Re- 
markably, the VBO is much larger than in the previous 
junction, although the change in the density between the 
constituent slabs is the same. The largest contribution 
to the VBO comes from the fact that in the hexagonal 
diamond the valence band width is larger than in cubic 
diamond, and the valence top has a splitting of 0.30 eV. 

The common feature of the semiconducting a-Si/c-Si 
interface and of crystalline Si/Si junctions is the reduced 
symmetry of one of the constituent slabs: this can be 
indeed considered as the origin of the VBO in these sys- 
tems, through a splitting of the valence band top edge 
manifold, more pronounced for the structure having the 
largest anisotropy (hexagonal diamond). 

Results for other a-Si/c-Si samples - We discuss now 
the results obtained for other a-Si/c-Si samples, which 
are characterized by a concentration of defects in the a-Si 
region at least of the order of ~10% and by a semimetallic 
behavior due to defect-induced gap states. Sample B 
is obtained following the same procedure for sample A, 
except that the starting configuration is obtained joining 
two separate cells (c-Si and a-Si, treated with the WWW 
model ||); sample C is obtained without the use of the 
WWW model, and following exactly all the steps (i) — (iv) 
previously listed, with a slow quench rate in step (i) in 
order to keep the number of defects very low; sample D 
is obtained with a higher quenching rate in step (i) and 
a TB finite-temperature annealing in step (iv); sample 
E with a TB athermic annealing; finally, sample F with 
with the use of EDIP-MD only (i.e., steps (Hi) — (iv) 
were skipped). 

Our key-result is that for all the samples from B to F 
the band alignment has almost the same value </> p =0.22 
eV, independent on structural details such as concentra- 
tion and distribution of defects and density of the amor- 
phous slab. The average electrostatic potential in the 
crystalline slab (V) c is very sensitive to the defect con- 
centration, and infact the values for the samples with 
intermediate defect concentrations (B, C) are very differ- 



ent from those with high defect concentrations (D, E, F), 
even by 0.48 eV (see Tab. I). Nonetheless, the Schottky 
barrier <f> p varies within a very small range (of 0.04 eV) 
for all these samples. 

This result for defect-rich samples is consistent with 
the general trend reported for the semiconductor/metal 
junctions pjj , i.e. that the Schottky barrier is poorly 
sensitive to the metal type but depends on the semicon- 
ductor. In fact, what changes from sample B to F is only 
the highly defected a-Si slab, which acts as the "metal- 
lic" constituent. Moreover, since c-Si is homopolar, 
neither interface-specific details could change the band 
alignment, as instead it may occurs in case of heterova- 
lent semiconductor heterostructures and metal/polar- 
semiconductor junction. The independence on interface 
details is also numerically verified from our calculations: 
in all our samples electric fields within the c-Si slabs are 
absent, whereas they would be present if the two inter- 
faces within the supercell were non-equivalent. 

Conclusions. - Present results 

for amorphous/crystalline Si interfaces can be rational- 
ized along the findings of the linear-response theory of 
the band alignment at semiconductor interfaces pl| and 
summarized as follows: (i) the band alignment is almost 
vanishing for defect-free samples and sizeable in the case 
of defect-rich samples; (ii) in the latter case it is inde- 
pendent on the particular defect concentration and distri- 
bution: such result can be ascribed to the semi-metallic 
behavior of the amorphous region, and explained by anal- 
ogy with the case of Schottky barriers which are almost 
independent on the metal type; (Hi) in any case it is not 
an interface-specific effect, but rather it is a bulk-related 
effect; (iv) by inspection of the different a-Si/c-Si samples 
and for comparison with other Si/Si crystalline junctions, 
the largest contribution to the band alignment is not due 
to changes in the density but rather to a broadening of 
the valence band and a splitting of its top manifold which 
typically occurs in the phases of reduced symmetry. 

This work was supported by the INFM Parallel Com- 
puting Initiative. 
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TABLE I. Main structural and electronic properties of 
the six different a-Si/c-Si samples considered in this work. 
The percentages of atoms in the amorphous part which are 
under-coordinated (Tjvjv<4) and over-coordinated (Tjvjv>4), 
where NN means "nearest neighbors", are indicated in the 
second and third columns. 



Sample 


TnN<4 


TnN>4 


c/a 


E F (eV) 


(V)c (eV) 


<fe (eV) 


A 


1% 


0% 


5.039 


6.26 


-0.11 


0.06 


B 


7% 


6% 


5.126 


6.22 


-0.30 


0.20 


C 


0.5% 


9% 


5.100 


6.33 


-0.24 


0.24 


D 


2% 


36% 


4.922 


6.71 


0.18 


0.21 


E 


3% 


27% 


4.922 


6.70 


0.15 


0.23 


F 


0.5% 


36% 


4.922 


6.68 


0.14 


0.22 




z(a ) 

FIG. 1. Profile of the electron density (dotted line, left 
scale) and of the macroscopic average of the total electro- 
static potential (solid line, right scale) for the sample D (see 
text). The position along the growth direction z is indicated 
in units of ao, the lattice parameter of the crystalline phase. 




FIG. 2. Scheme of the band alignment at a-Si/c-Si interface 
with defect-rich a-Si, with indication of: a typical profile of 
the electrostatic potential in the supercell (thick solid line), its 
average over the whole supercell ((V)tot), its average in the 
crystalline region ({V) c ), the Fermi energy of the supercell 
(Ef), the valence band top edge of c-Si (E VtC ) with respect 
to (V)tot, and finally the Schottky barrier <j> p . 
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